Klein tunneling in carbon nanostructures: a free particle dynamics in disguise 
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The absence of backscattering in metallic nanotubes as well as perfect Klein tunneling in potential 
barriers in graphene are the prominent electronic characteristics of carbon nanostructures. We show 
that the phenomena can be explained by a peculiar supersymmetry generated by a first order 
Hamiltonian and zero-order supercharge operators. Like the supersymmetry associated with second 
order reflectionless finite-gap systems, it relates here the low-energy behavior of the charge carriers 
with the free-particle dynamics. 
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Graphene is a genuine two dimensional material com- 
posed of the carbon atoms that form a honeycomb lattice. 
Three of the valence electrons of each carbon atom partic- 
ipate in the interatomic interaction, while the fourth one 
contributes to the conductivity of the crystal. Graphene 
has been studied theoretically for a long time, see e.g., 
[H, Q. However, its experimental observation trig- 
gered a real boom of both theoretical and experimental 
analysis 0, @, 0. 

The material manifests extraordinary electronic prop- 
erties, which are the consequence of an unusual dynamics 
of the low-energy charge carriers. It was pointed out in 
Q that the tight-binding description of the system is re- 
duced to the massless Dirac equation in the low-energy 
approximation. This makes graphene an ideal test field 
for (2 + l)-dimensional QED [fj|; due to the low Fermi 
velocity vf, cjvp ~ 300, it is possible to simulate rela- 
tivistic effects in condensed matter systems which would 
be unreachable experimentally otherwise. 

It was predicted , that the scattering of the rela- 
tivistic electrons on the potential barrier is qualitatively 
different from the nonrelativistic case. The particles can 
tunnel the barrier without reflection, provided that its 
height tends to infinity. This is in contrast to the nonrel- 
ativistic regime where the tunneling would be exponen- 
tially suppressed [§[ . This phenomenon, known as Klein 
tunneling, is not experimentally realizable with elemen- 
tary particles nowadays due to the extreme electric field 
needed to observe the predicted difference between rela- 
tivistic and nonrelativistic scattering Q. 

The scattering of the low-energy quasiparticles in 
graphene on the barrier with translational symmetry in 
one dimension was analyzed in 0, 0, HH- The absence 
of backscattering was noticed for normal incidence. The 
effect is independent of the height of the barrier and, 
hence, is testable experimentally (l2| . A similar phe- 
nomenon was observed earlier [13[ and discussed theoret- 



despite the presence of a scattering potential generated 
by impurities. The absence of backscattering was under- 
stood as a consequence of topological singularity identi- 
fied with a Dirac point, see 17 1, [18J, or as a result of the 
pseudospin conservation 0. 

We provide here a simple, alternative explanation for 
the absence of backscattering in the carbon nanostruc- 
tures within the framework of super symmetric quantum 
mechanics. We shall discuss a broad class of potentials 
in graphene as well as in the metallic nanotubes with the 
range exceeding the interatomic distance. 

The honeycomb lattice is a superposition of two tri- 
angular sublattices, A and B. The eigenstate ^ of the 
Hamiltonian can be then written as * = ca^ a + cb^ b , 
where ^a and b are atomic wave functions of the sub- 
lattices whereas ca and cb are slowly varying amplitudes. 
^ is a Bloch function, which acquires a nontrivial phase 
factor when shifted by a translation vector R of the Bra- 
vais lattice, \l/(k, x + R) = e tkR \l/(k, x). Fermi surface of 
graphene is formed by discrete points. There are six of 
them in the first Brillouin zone, situated in its corners, 
see FiglTJ In the analysis of the low-energy behavior of 
the charge carriers, it is sufficient to consider just two of 
them, denoted as Dirac points K and K' = K. The re- 
maining four Dirac points do not represent distinct elec- 
tronic states. They can be obtained either from K or K' 
by translations in the reciprocal lattice. 

In the vicinity of Dirac points, the behavior of the sys- 
tem is described by the massless Dirac equation. When 
the effective Hamiltonian is considered in the valley of 
the point K and expanded up to the terms linear in the 
momenta, the energy eigenvalue equation acquires the 
form (we put h = 1) 



Hip = -i v F (crid x + a 2 d y ) ip = Eip, 



(1) 



ically [14| , [15| , [16] in the context of electron transport 



in carbon nanotubes. The perfect transmission of the 
low-energy charge carriers occurs in metallic nanotubes 



where o~i t 2 are Pauli matrices. Spinor -0 reads explicitly 
tp = (ca, cs)'e I,5kx , where 5k = k — K and t is a trans- 
position. Degree of freedom associated with the compo- 
nents ca and cb is called pseudospin, see 0. 

A single wall carbon nanotube can be created by rolling 
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Figure 1: (color online). The hexagon represents 
the first Brillouin zone of graphene. The blue con- 
tours show constant values of the energy modulus \E\, 
E — ±7oy3 + 2coskai + 2coska2 + 2cosk(ai — a2), 70 = 
2(V3a)~ 1 VF- The primitive translation vectors of the hexag- 
onal lattice are fixed as ai = a(0, 1) and a2 = a/2 (\/3, l), 
while K = a _1 (0, 47r/3), a = V%d, d is the nearest bond 
length of the graphene lattice (we put here a = 1 and 
vf = V3/2). The deep blue corresponds to the lowest values 
of \E\. The energy surface forms deep valleys near the Dirac 
points. The parallel lines illustrate the allowed momenta for 
the zig-zag nanotube (6,0), i.e. kCh/27r £ Z, 



where the value of e S |o, ± 3 | ^ h | j depends on the type 
of the nanotube : it is zero for metallic nanotubes and 
nonvanishing for semiconducting ones. The parameter e 
can be alternatively regarded as a fictitious magnetic flux 
22] or the mass of the quasiparticle. Thus, the Hq, given 
by @ with e = 0, coincides with the Dirac Hamiltonian 
of the free massless particle in one dimension. 

The real- world nanostructures are not perfect. The 
crystal can have impurities; or it can be in presence of 
an external field. We suppose that these effects are rep- 
resented by a potential V, which is vanishing at infinity, 
varies smoothly on the scale of the interatomic distance 
and is of the range larger than it (and, hence, does not 
cause the intervalley transition of the electronic states). 
Then it can be consistently incorporated into (p} and ([2]) 



as the unit matrix multiplied by V 15]. The station- 
ary evolution equation for metallic nanotubes acquires a 
form 



Hyip = (-iv F aid x + V(x))ip — Eip. 



(3) 



The same Hamiltonian Hy can be obtained as a reduc- 
tion of the Hamiltonian of graphene with a potential bar- 
rier, H + V, as long as potential is translationally invari- 
ant in the y-direction, i.e. V = V(x), 



H v = e 



-iSk y y 



(H + V)e 



iSk y y 



5fc„=0- 



(4) 



up and gluing appropriately a single graphene strip [l9| . 
The circumference (chiral) vector Ch is an integer multi- 
ple of the primitive translation vectors ai and a.2 of the 
Bravais lattice, Ch = niSLi + n^a^.- It defines uniquely 
the nanotube [20 ] by the periodic condition imposed on 
the wave functions, ^P(x + Ch) = ^(x). Let us sup- 
pose that the coordinate system is chosen such that 
the y axis points in the direction of the chiral vector, 
k Ch = k y | Ch | • Taking into account Bloch character of 
the wave functions, one can see that ^ has the required 
periodicity as long as k Ch = k y \ Ch \ = 27rg for an integer 
q. Hence, momentum k is quantized in the y-direction 
and the allowed values form equidistant lines in the first 
Brillouin zone, with a step jg^j, see FigQ]. 

There are two qualitatively different positions of the 
Dirac point K with respect to the allowed lines, in co- 
herence with two main classes of carbon nanotubes. In 
case of metallic nanotubes, K is crossed by one of the 
lines. For semiconducting nanotubes there is a minimum 
distance mm|k — K = 27r/|3Ch| between the Dirac point 



and the closest line [21[, which leads to the opening of 
the gap between the valence and conduction bands. 

The low-energy behavior of the charge carriers in the 
nanotube with the chiral vector Ch is approximated by 

Help = vf (-i(?id x + ea 2 ) ip = Eip, (2) 



In this case, Hy represents the energy operator of the 
particle with normal incidence (dk y = 0) on the barrier. 
The Hamiltonian Hy describes reflectionless system, in- 
dependently on the particular form of the potential V. 
This peculiar fact was explained by the presence of Berry 
phase in the wave functions, which causes the Born series 
for backscattering to vanish identically 



17] 



Let us look, however, at the problem from a new per- 
spective. In nonrelativistic quantum mechanics, the re- 
flectionless systems play an important role in the theory 
of solitons and are intimately related to the nonlinear in- 
tegrable systems (23[. Their particular properties come 
hand in hand with the existence of the supersymmetry, 
that is based on the Darboux-Crum transformations [23| . 
241 ] , 25j . The supercharges intertwine such systems with 
the free-particle model and stay behind the absence of 
backscattering in the conduction band. 

This rises the question whether the perfect tunneling in 
the carbon nanostructures has a similar algebraic back- 
ground. The answer is affirmative, despite the fact that 
the relativistic Hamiltonian Hy is of the first order, con- 
trary to the second order Hamiltonians of the nonrela- 
tivistic reflectionless systems [26j |. 

Let us define the hermitian operators 



H = 



Hy 
H 



r = 



1 

-1 



(•5) 
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where 1 is the unit two-by-two matrix. Besides the op- 
erator r, the Hamiltonian % has two other symmetries, 



o W 

U 



U 2 = iTU 1 , 



(6) 



where U = U (x) is a unitary operator of a local chiral 
rotation (chiral gauge transformation), 



l sin a o\ , 



(7) 



dependent on the interaction potential, 

a(x) = — ( V(t)cLt. (8) 

VF J 

They satisfy relations 

[U,U a ]=0, {U a ,Ub} = 2S ab l, a, b = 1,2. (9) 

The anticommutator of IA\ and U 2 is a zero-order poly- 
nomial in %, to be proportional to the central element 
1 which is the unit 4x4 matrix. The grading opera- 
tor T classifies the operators dU and 1 as bosonic (they 
commute with T) while operators are fermionic (they 
anticommute with T). 

Relations ([9]) constitute the N = 2 zero-order super- 
symmetry extended by the central charge 1 and graded 
by r. Like in the nonrelativistic case of a reflectionless 
system with the n-gap, second order Hamiltonian (where 
the order n of supersymmetry is fixed by the number n of 
bound states [25|), this structure underlies the absence 
of the backward scattering in the system given by the 
first order Hamiltonian Hy. The relation UHy = HqU, 
implied by the commutator in ([9]) and the unitarity of U, 
reveal the unitary equivalence of Hy with the free mass- 
less Dirac Hamiltonian Ho- This proves the absence of 
backscattering : the setting given by Hy is unitary equiv- 
alent to the free massless particle system and, hence, it 
shares its trivial scattering properties. 

Because of this peculiar supersymmetric structure, 
based on the nontrivial unitary equivalence, all the inte- 
grals of the free-particle system represented by H have 
their analogs in the system given by Hy. In particu- 
lar, the pseudospin is conserved in Hy just because of 
[Ho,ai] = [U,(Ti] = 0. The free-particle momentum 
—id x transforms into the integral — id x + -^<j\V{x) = 
—iD x for Hy, that has a form of a covariant deriva- 
tive with a chirality-dependent charge. Then Hy, like 
Ho, is presented as a composition of the integrals, Hy — 
-ivFO\V x . 

The unitary equivalence of Hy with the free-particle 
system is broken when the nonvanishing effective mass 
(the coefficient of cr 2 ) is present in the impurity Hamilto- 
nian. It emerges in semiconducting nanotubes {Hy = 
H e + V with m — e — ±27r/|3Ch|), or when 
other than normal incidence of the particles is con- 



V)e m w\ Skv=rnji0 = H m + V, where H m = H + v F ma 2 . 
In these cases, the unitary transformation of Hy yields 
UHyU^ 1 = H + mv F cr 2 (x), a 2 (x) = Uo 2 U~ x = 
cos 2a a 2 — sin 2a 03 , and 



UHy = H m U — 2v F msuia a^. 



(10) 



Therefore, the scale of supersymmetry breaking in the 
massive case is of the order of m, and the contribution 
of the potential is controlled by the factor |sina < 1. 
Hence, for the close-to-the-normal incidence (m = 6k y ~ 
0), the potential barrier remains almost perfectly trans- 
parent for any V(x). This is coherent with [Toj . where 
the scattering on n—p junction was analyzed. For general 
values of the effective mass, the scattering properties of 
Hy are nontrivial, however, and depend on the explicit 
form of the potential. The quasiparticles can be confined 
in graphene, the charge carriers get localized in semicon- 



27]. 



sidered in graphene system (U), H\ 



-i8k y y 



ducting nanotubes, see [16 

Zero-order supersymmetry, based on a unitary equi- 
valence of the superpartner Hamiltonians, can be for- 
mally constructed for any quantum system. It is suf- 
ficient to make a unitary transformation of an initial 
Hamiltonian to get its superpartner. However, the su- 
persymmetry is nontrivial and manifests its predictive 
power when the Hamiltonians describe different physics, 
like in ([5]). There, it provided a simple explanation for the 
absence of backscattering in the considered carbon nanos- 
tructures. As an example of another application, it may 
be used to determine the s-wave of the Dirac operator in 
polar coordinates [f|, H(r,9) = H (r) — ivFcr 2 ^dg + V(r), 
where H (r) — ivp(Ti(d r +'i-), for any potential V(r). In- 
deed, H(r, 9) acts on the subspace of s- waves (ip = ip(r)) 
as Ho(r) + V(r). We can use an r-dependent unitary 
mapping to get rid of the interaction term and to obtain 
the operator Ho(r) with analytically computable eigen- 
states. 

In the context of graphene, the supersymmetry ap- 
peared earlier in the analysis of the zero-energy states in 
presence of external magnetic field [28[ , and in the study 
of the quantum Hall effect in particular [29[ . In that case, 
the first order supercharges are proportional to the Dirac 
Hamiltonian, while the supersymmetric Hamiltonian is of 
the second order. The structure in ^ is completely dif- 
ferent: the supersymmetric Hamiltonian T-L is of the first 
order and the existence of the zero-order supercharges 
lA a provides a complete information on the eigenstates 
of Hy. Together with (fTTJ| . it suggests to analyze the 
system with small effective mass perturbatively. 

The algebraic framework presented here can be broad- 
ened in different ways. For instance, the apparent simi- 
larity with the nonrelativistic reflectionless systems sug- 
gests to extend the analysis by employing the Darboux- 
Crum transformations in the context of the higher-order 
(nonlinear) supersymmetry (2(| |3(| . We believe that in 
such a generalized form the supersymmetry can serve as 
a useful tool in the study of the low-energy excitations 
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of the charge carriers in the wrapped graphene, where 
the dynamics is governed by the Dirac Hamiltonian in a 
curved space. Further discussion on this problem, and 
other possible applications indicated above goes beyond 
the scope of the present article. 
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